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Strong magnetic-ferroelectric coupling in hexagonal HoMnOs has been observed 
previously at spin reorientation temperature between P6 3 cm and the P6 3 c m mag- 
netic phases. In contrast, YM11O3 has only a single magnetic phase (P6 3 cm') and no 
sign of strong magnetic-ferroelectric coupling. In order to investigate the P6 3 cm' to 
P6 3 c'm spin reorientation transition, single crystals of Hoi_ x Y x Mn03 at varying 
compositions were grown. Neutron scattering measurements reveal quasielastic scat- 
tering in the YMnC>3-like phase (P6 3 c'm) centered on the Bragg peak correspond- 
ing to HoMn03-like phase (P6 3 cm') and vice versa. This scattering may be due 
to short-lived fluctuations into the "wrong" magnetic phase at domain boundaries. 
These results suggest that there is strong phase competition between the P6 3 cm' and 
P6 3 c'm phases even in pure YMn03, and this competition may play an important 
role in the magnetic- ferroelectric coupling observed in HoMnOs. 
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I. INTRODUCTION 

The hexagonal manganite multiferroics (i?Mn0 3 , R = Rare earth) are a family of isostruc- 
tural compounds which order ferroelectrically at high temperatures and antiferromagnet- 
ically at low temperatures [1]. Magnetic order in this family is created by 2 layers of 
Mn 3+ stacked on top of each other at Z=0 and Z=l/2 with in-plane 120° angle between 
neighboring Mn 3+ S — 2 magnetic moments within each plane [2]. Much of the inter- 
est in this family of multiferroics is due to the strong magnetic-ferroelectric coupling ob- 
served in HoMn0 3 [3-6]. As HoMn0 3 undergoes a spin reorientation transition between a 
P6 3 c'm magnetic phase and a P6 3 cm' phase with decreasing temperature [7], there is a very 
large increase in the c-axis dielectric constant at the spin reorientation transition temper- 
ature Tsr [3, 4] along with a change in the magnitude of the ferroelectric polarization [5]. 
HoMnC>3 is the only pure member of this family of multiferroics where this strong magnetic- 
ferroelectric coupling at a spin reorientation transition has been observed, and the physics 
behind this strong coupling is not well understood. In contrast, YMn0 3 orders antiferro- 
magnetically in P6 3 cm' phase symmetry and does not undergo any spin reorientation at zero 
field [1, 7]. 

The different magnetic phases in i?Mn0 3 are distinguished by how spins in neighboring 
planes are correlated, and the determining factor seems to be the position of the Mn 3+ 
ion within the unit cell respect to a critical value of 1/3 [8]. The mechanism responsible 
for the large magnetic-ferroelectric coupling observed in HoMn0 3 is still not well under- 
stood, though. One approach to try to better understand this coupling is to study how 
various parameters change the spin reorientation temperature where this coupling is seen. 
Such work has included magnetic fields [1, 3, 4], electric field [6], pressure [9], and chemical 
substitution [10-12]. By doping either Er or Y into HoMn0 3 , the spin reorientation transi- 
tion temperature Tsr can be decreased or increased, respectively. This makes sense, since 
ErMn0 3 orders in the P6 3 c'm phase and YMn0 3 orders in the P6 3 cm' phase. However, we 
have previously found that the phase boundary between these two phases as a function of 
rare earth concentration is actually complex, with qualitatively different behavior for other 
Ri_~Y-Mn0 3 compounds [13, 14]. 

Previous work done by H.D Zhou et al. investigated the phase diagram of Hox_2.Yj.MnO3 
with magnetic susceptibility, specific heat and dielectric constant [11]. For x > 0.9, 
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Hox_j.Y2.MnO3 orders in the P6 3 cm' phase and no spin reorientation occurs. They at- 
tribute the magnetic susceptibility and specific heat increases at Tsr to partial ordering of 
Ho 3+ moments, both of which are suppressed upon Y doping. The dielectric susceptibility 
anomaly at Tsr was observed for all samples with a spin reorientation transition, even up 
to x — 0.8. Pure YMn03 does not have any dielectric susceptibility anomalies along the 
c axis [15], not even the slight kink observed at the Neel temperature Tjy in HoMn0 3 [3]. 
In order to further study the evolution of magnetic order in Hoi_, z Y :r Mn03, we have per- 
formed neutron scattering measurements of the spin dynamics in single- crystal samples. 
Our results may help illuminate the mechanisms responsible for the coupling of magnetic 
and ferroelectric order in HoMn03. 



II. EXPERIMENT 



Single crystal samples of Hoi^^Y^MnOs were grown using an optical floating zone tech- 
nique, where x — 0, 0.65, 0.75, 0.825, 0.88, 1. Powders of Ho 2 3 , Y 2 3 and Mn0 2 were mixed 
and ground together using a pestle and a mortar and calcinated at 1200 °C for 8 hours. We 
repeated this process thrice to make it uniform. Afterwards we reground and packed in 
to both seed and feed rods and sintered at 1450 °C for 8 hours. Crystals were grown in 
a 4-mirror floating zone furnace under air at atmospheric pressure. Grown crystals were 
approximately 2 to 5 grams. Neutrons scattering of single crystal samples were performed 
at the Missouri University Research Reactor (MURR) using the triple-axis spectrometer 
TRIAX, using a closed-cycle helium displex to cool the samples to low temperature. 



A. Elastic Neutron Scattering 

When the spins undergo 90° reorientation from P6 3 c'm to P6 3 cm' there is a change in 
intensity of the magnetic Bragg reflections. In the P6 3 c m phase between T^v and Tsr the 
(1,0,1) Bragg reflection dominates, while in the P6 3 cm' phase below Tsr and in samples with 
no spin reorientation the (1,0,0) Bragg reflection dominates [16]. The (1,0,2) reflection is 
both structural and magnetic, and has significant scattering but different intensities in both 
magnetic phases. Figure 1 shows example order parameter scans for Hoo.75Y . 25 Mn0 3 . The 
Bragg scattering intensity scales with the square of the antiferromagnetic order parameter 



5 



M S t- The Neel temperature TV was obtained by fitting the peak intensity / of the dominant 
reflection above 60 K to the form of a second-order phase transition: 

/ oc Ml oc (T N - T) 2/3 . (1) 

The fits include a small Gaussian spread in the transition temperature (shown by the error 
bars for Tjy )• The critical exponent j3 value varied from 0.187 to 0.333 and was very sensitive 
to the fitting range of data. These critical exponent values should not be considered as the 
true critical exponents of the Neel transition, but T^ obtained from fitting can be considered 
reliable since they were not very sensitive to the fitting range. The shift in scattering intensity 
from the (1,0,1) to the (1,0,0) with decreasing temperature indicates Tsr, and error bars 
indicate the width of this first-order transition. With increasing Yttrium concentration, 
T S r increases towards to T^ and eventually at x — 0.88 T S r completely vanishes. Results 
for T N and Tsr are shown in Figure 1. The magnetic phase diagram for Hoi-^Y^MnOs we 
obtained is consistent with previous work done by Zhou et al. [11]. 



B. Inelastic Neutron Scattering 

Below the Neel temperature, inelastic measurements of YMnC>3 revealed quasielastic 
scattering mainly around (1,0,1), and secondary scattering (1,0,2) Bragg positions. Example 
data taken at 30K are shown in Figure 2, which revealed a peak centered at zero energy 
transfer, but with an energy width much larger than instrumental resolution. Scans along the 
H and L directions (in-plane and out-of-plane, respectively) were performed at a constant 
energy transfer of 0.7 meV in order to reduce elastic Bragg scattering. Previous work done by 
Sato et al. [17] reported quasielastic scattering with similar in-plane H momentum transfers, 
but with no observed out-of-plane momentum dependence. 

The data were fit to a structure factor S(Q,u) convoluted with the spectrometer resolu- 
tion, where the structure factor model was given by the following equation: 

1 Tu Wf Wf ^4 

W\\ is the in-plane width of the scattering, W± is the out-of-plane width, and T is the energy 
width. AH, AK, and AL are the difference in H, K, and L respectively from the center of the 
quasielastic scattering. The first fraction is the Boltzmann term, followed by Lorentzians 
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in energy and in-plane momenta H and K. The last term is a Gaussian in out-of-plane 
momentum L. Figure 2(c) shows fits using both the above static structure factor (solid line) 
and a similar form using a Lorentzian instead of a Gaussian for the out-of-plane dependence, 
and the Gaussian fits produced visibly better results. The out-of-plane correlation length in 
lattice units is given by 2/W±, and results are plotted in Figure 3. The data were fit to the 
heuristic form 

Ik(T n -T) 2X . (3) 

where A 0.26 for the above fit. With increasing temperature the out-of-plane correlation 
length decreases, and above T^r no apparent out-of-plane correlations persist. The total 
quasielastic scattering intensity also decreases with decreasing temperature below T/v- 

Similar quasielastic scattering was also observed in Hoi-^Y^MnOs samples for < x < 
1. Figure 4 shows scattering intensity maps as a function of H and L at three different 
temperatures for Hoo.25Yo.7sMn03. Measurements were taken at a constant energy transfer 
of 0.5 meV to reduce elastic background scattering. Above T^v, the scattering was broad in 
the H direction and had no L dependence, as shown in Figure 4(a), indicating short-range 
in-plane correlations with no out-of-plane correlations. Near Tjv, the quasielastic scattering 
showed two quasielstic peaks (Figure 4(b)) centered at the (1,0,0) and (1,0,1) positions. 
Below T^, scattering became sharper in both in plane and out of plane directions, and was 
only centered around the (1,0,0) position at 60K (Figure 4(c)). 

For Roi- x Y x Mia0 3 samples, which undergo a spin reorientation transition, the location 
of the quasielastic peak was observed to shift at this transition. Figure 5 illustrates this 
shift from (1,0,0) to (1,0,1) at T S r fc> r two different concentrations of Hoi^Y^MnOs. When 
the system is in the P6 3 cm' phase below Tsr, the magnetic Bragg reflection is located 
at (1,0,0) (lower panel) and the quasielastic scattering is located at (l,0,l)(upper panel). 
Above Tsr the magnetic Bragg reflection is located at (1,0,1) and the quasielastic scattering 
shifts to the (1,0,0) position. Quasielastic scattering was also observable at the (1,0,2) 
position in both phases. These results match similar observations in Erx_a.Yj.MnO3 an d 
D yi __Y_Mn0 3 [13, 14]. 
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III. DISCUSSION 



Diffuse neutron scattering in i?MnC>3 similar to our observations has been previously 
reported in several papers [17-21]. Sato et al. reported quasielastic scattering above and 
below Tjy in pure YMn03, with an in-plane width that became resolution limited below 
TV, an energy width that remained wider than resolution down to the lowest temperatures, 
but with no reported out-of-plane correlation [17]. Park et al. also reported quasielastic 
scattering above and below in YMnC>3, but their samples were powder and so the data 
does not provide directional information [18]. Roessli et al. measured quasielastic scattering 
in pure YMn0 3 around the (1,0,0) position above T N but not below [19]. Far above Tjy, the 
scattering had no out-of-plane correlations, but near TV the scattering developed L depen- 
dence with a peak centered around (1,0,0), and a peak shape well-described by a Lorentzian 
in both the in-plane and out-of-plane direction. Polarized neutron measurements indicated 
that the quasielastic scattering was due to in-plane fluctuations. Lewtas et al. reported 
diffuse elastic scattering in LuMn0 3 around the (1,0,0) position above TV, while the pri- 
mary magnetic Bragg reflection below T N was at the (1,0,1) position [21]. This scattering 
may have been quasielastic, but no energy-dependent measurements on the diffuse scatter- 
ing was reported. Sekhar et al. reported diffuse scattering below TV in powder neutron 
measurements of Eri_, E Y. r Mn03, but only in samples with high Yttrium concentrations [20]. 
This scattering occurs at 29 angles similar to the (1,0,0) and (1,0,1) peaks, but they did not 
identify which of those two positions this scattering was associated with. 

There is agreement among all the single-crystal measurements that well above T/v, the 
diffuse/quasielastic scattering has no out-of-plane correlations and very short in-plane cor- 
relations. Both Lewtas et al. and Roessli et al. report that the scattering develops L 
dependence centered around the (1,0,0) position in LuMn0 3 and YMn0 3 , respectively, but 
neither report measuring around the (1,0,1) [19, 21]. This omission is significant because the 
(1,0,1) position is the magnetic Bragg position for LuMn03, but the (1,0,0) is the magnetic 
Bragg position for YMnC>3. Below the Neel temperature, there is significant discrepency be- 
tween our work and the work of Sato et al. [17], since the latter report that the quasielastic 
scattering has no out-of-plane dependence, while we observe very strong out-of-plane de- 
pendence below T N . We speculate that the Sato et al. measurements may simply not have 
extended far enough. Only one L scan was shown, and the range covered, from L = —0.2 to 
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L = 0.5, corresponds to an essentially flat region of our scans, as seen in Figure 2. Our own 
work on Hoi_.jY,2.Mn0 3 , as well as our previous work on Er^Y^MnOs and Dy 1 _ x Y a ,Mn0 3 , 
has shown that there are significant out-of-plane correlations in the quasielastic scattering 
of i?Mn0 3 below TV- However, while these quasielastic excitations have now been observed 
in YMn0 3 [17-19] and LuMn0 3 [21], as well as Hoi-^Y-MnOa, E_i__Y_Mn0 3 [13, 20], 
and Dyi_j.Yj.MnO3 [14] at high Yttrium concentrations, they were not observed in inelas- 
tic neutron scattering measurements of pure HoMn0 3 [22] or in pure ErMn0 3 diffraction 
measurments [20]. 

Well above T w , the quasielastic scattering in these materials shows no sign of L depen- 
dence. Correlations are in-plane only, and so no distinction can be made between phases 
which differ in their out-of-plane stacking arrangement. Near Tjy, however, L dependence 
begins to emerge as correlations begin to develop between MnO planes. The presence of 
quasielastic scattering at both the (1,0,0) and (1,0,1) positions near the T^ indicates that 
correlated regions of both phases exist above T/v- Neither phase is strongly preferred, but 
instead they appear to compete with each other. According to our interpretation there 
should be a peak in the quasielastic scattering just above Tjy at the (1,0,1) position in 
LuMn0 3 in addition to the observed peak around (1,0,0) [21], but it may simply not have 
been measured. 

Below T N , YMn0 3 orders in only one phase, and quasielastic scattering for the or- 
dered phase disappears as the correlations become static. But quasielastic scattering for 
the "wrong" phase persist below T^v- These patches are broad in-plane and narrow out-of- 
plane. We hypothesize that the quasielastic scattering originates from domain boundaries 
within our crystal. Mn 3+ moments located between different domains will be frustrated. 
At c-axis domain boundaries, the weak, frustrated out-of-plane interactions may prevent 
long-range static order, but in-plane correlations within those boundaries could still lead to 
short-range correlated fluctuations. Thermal fluctuations may excite spins at the boundary 
into the wrong magnetic phase, which still satisfies the in-plane interactions. If the cor- 
relations which form are short-lived, they will lead to quasielastic scattering. Because the 
in-plane coupling is so much stronger than the out-of-plane coupling, fluctuations at in-plane 
domain boundaries may not be able to form correlated patches of any significant size, and 
so would not contribute to the observed quasielastic scattering. This would also explain 
the line shape of the quasielastic scattering: the shape of the domain walls would constrain 
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correlations along the c axis, but correlations could grow to much larger size along the in- 
plane direction. The neighboring ordered domains will not only act to limit the size of the 
out-of-plane correlation length, they may also cause those correlations to decay faster than 
an exponential. A Lorentzian line shape indicates correlations which decay exponentially 
with distance, while a Gaussian indicates correlations which decay much faster than expo- 
nential. This contrasts with the analysis of Roessli et al, who used a Lorentzian to analyze 
the L dependence of the quasielastic scattering above T^r [19]. Without an ordered phase 
to constrain the fluctuations, correlations could extend as a decaying exponential along the 
c axis, so their model is still consistent with our hypothesis. 

We also hypothesize that the existence of quasielastic scattering in Yttrium-rich i?Mn0 3 but 
not pure HoMn0 3 is due to the larger in-plane easy- axis anisotropy in HoMn0 3 . This in- 
plane anisotropy for Mn 3+ moments gives rise to a spin gap which was measured by Fabreges 
et al. for both Yttrium and Holmium, who found a significantly larger gap for HoMn0 3 than 
for YMn0 3 [8]. This anisotropy acts as an energy barrier to in-plane rotation, and so a 
large anisotropy may prevent the quasielastic fluctuations at the domain boundaries. Fur- 
ther evidence for this interpretation comes from the polarized scattering measurements of 
Roessli et al, who found that the quasielastic scattering is due to spins fluctuating within 
the plane [19]. Out-of-plane fluctuations gave rise to inelastic scattering due to the easy- 
plane anisotropy, which is much larger than the in-plane anisotropy for all i?Mn0 3 . The 
observation of quasielastic scattering in LuMn0 3 is consistent with this interpretation, since 
the spin-wave dispersion measured by Lewtas et al. had a vanishingly small gap indicating 
very little in-plane anisotropy [21]. 

One possible objection to our interpretation is that the observed quasielastic scattering 
may simply be the spin-wave dispersion approaching zero energy. The spin-wave measure- 
ments by Fabreges et al. [8] do indicate that the dispersion is close to gapless for YMn0 3 at 
the (1,0,1) position, and in HoMn0 3 it has a minimum at either the (1,0,1) or the (1,0,0), 
depending on the temperature, in the same place as the observed quasielastic scattering 
for YMn0 3 and Hoi^Y^MnOs. However, the observation of quasielastic scattering at the 
(1,0,2), where the dispersion should move to higher energy in YMn0 3 due to out-of-plane 
coupling, indicates that the quasielastic scattering is not simply a measurement artifact of 
a gapless dispersion, but is instead related to correlations matching an ordered magnetic 
phase, albeit the wrong phase. 
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Our model of quasielastic scattering in i?Mn0 3 explains both our own observations as 
well as a number of disparate previous results. It also explains the absence of quasielastic 
scattering for HoMn0 3 . If our model is correct, then other RMn0 3 compounds with large 
in-plane anisotropy, such as YbMn0 3 [8], should not exhibit quasielastic scattering either. 
This scattering indicates strong competition between competing 3D-ordered phases. It is 
particularly interesting that even though HoMn0 3 exhibits both P6 3 cm' and P6 3 c'm phases, 
quasielastic scattering is absent. The in-plane anisotropy appears to stabilize both magnetic 
phases in HoMn0 3 . This suggests that the spin reorientation transition in HoMn0 3 is driven 
by something besides purely magnetic interactions. We know that the magnetism and the 
ferroelectricity are coupled at this transition, and a better understanding of the dynamics 
involved may help illuminate the nature of this magnetic-ferroelectric interaction. 
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FIG. 1. Order parameters for H00.75Y0.25M11O3 concentration. In the left panel, white squares 
show the (1,0,0) peak intensity, gray triangles show the (1,0,1) peak intensity, and black circles 
shows the intensity of (1,0,2) peak. The (1,0,0) reflection is associated with the PG^c'm phase, 
the (1,0,1) reflection is associated with the P6 3 cm phase, and the (1,0,2) relfection is associated 
with both phases but at different intensities. Lines show fits to the (1,0,1) order parameter data 
near TV. The right panel show the phase diagram extracted from order parameter measurements. 
White circles indicate Ttv, black triangles indicate Tsr- 
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FIG. 2. Measurements of the quasielastic scattering in YMnC>3 at 40K. Scattering is centered at 
the Q=(1,0,1) position, and scans were taken with varying energy (a), in-plane momentum (b), and 
out-of-plane momentum (c). Momentum scans were taken with fixed energy transfer of 0.7 meV to 
reduce elastic scattering background. Note that quasielastic scattering is also present around the 
(1,0,2) position. Lines are fits to the data described in the text. The solid line in (c) shows the fit 
result using a Gaussian for the out-of-plane dependence, while the dashed line shows the fit result 
using a Lorentzian. 




Temperature (K) 

FIG. 3. Out-of-plane correlation length for quasielastic scattering in YMnOs versus temperature. 
The line is a fit to the heuristic form described in the text. 
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FIG. 4. Ql v Qh scattering intensity maps for Hoo.25Yo.75Mn03 at a constant 0.5 meV energy- 
transfer. Scans were taken (a) at 80K (above TV), (b) 72K (approximately TV) and (c) 60K 
(below TV). Above T^r, the quasielastic scattering shows no out-of-plane dependence. At T^r, the 
quasielastic scattering begins to concentrate around both the (1,0,0) and (1,0,1) positions. At 60K, 
the scattering is centered around the (1,0,0) position, although the magnetic Bragg scattering is 
located at the (1,0,1) position at this temperature. 
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FIG. 5. Top panels: Measurements of the out-of-plane dependence of quasielastic scattering versus 
temperature for two compositions of Hoi_ x Y x Mn03. Bottom panels: Order parameter measure- 
ments of the same samples over the same temperature range. At the spin reorientation temperature, 
elastic Bragg scattering shifts from the (1,0,1) position to the (1,0,0) position upon cooling, while 
quasielastic scattering shifts from the (1,0,0) position to the (1,0,1) position at the same temper- 
ature. Scattering with very narrow L width is also visible in the top panel, and is due to the tail 
end of the resolution ellipsoid clipping the Bragg reflection 
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